Entropy current for the relativistic Kadanoff-Baym equation and H-theorem in 0{N) theory with 

NLO self-energy of expansion 
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We derive an expression of the kinetic entropy current in the nonequilibrium 0{N) scalar theory from the 
Schwinger-Dyson (Kadanoff-Baym) equation with the 1st order gradient expansion. We show that our kinetic 
entropy satisfies the H-theorem for the leading order of the gradient expansion with the next-to-leading order 
self-energy of the expansion in the symmetric phase, and that entropy production occurs as the Green's 
function evolves with an nonzero collision term. Entropy production stops at local thermal equilibrium where the 
collision term contribution vanishes and the maximal entropy state is realized. Next we also compare our entropy 
density with that in thermal equilibrium which is given from thermodynamic potential or equivalently 2 particle 
irreducible effective action. We find that our entropy density corresponds to that in thermal equilibrium with 
the next-to-leading order skeletons of the 1 /N expansion if skeletons with energy denominators in momentum 
integral can be regularized appropriately. We have a possibility that memory correction terms remain in entropy 
current if not regularized. 

PACS numbers: 



I. INTRODUCTION 

The 2 Particle Irreducible (2PI) effective functional tech- 
nique provides a powerful tool to deal with controlled 
nonequilibrium dynamics with non-secularity and late time 
universality |[ll01- It gives the Kadanoff-Baym (KB) equa- 
tion which describes the dynamics for quantum fluctuations of 
fields. In 1960s, Baym and Kadanoff studied the Schwinger- 
Dyson (SD) equation for the two-point function G{x, y) for 
nonequilibrium systems jslH] on the basis of a functional ap- 
proach developed by Luttinger and Ward |5i]. Then Baym 
reformulated it in terms of variational principle, introducing 
the so-called <I>-derivable approximation which is given 
by a truncated set of closed 2PI diagrams. The main virtue of 
this approximation is that the resulting equations conserve the 
charge, energy and momentum of the system. This approach 
has been reformulated by use of path integral method IIT HlOtl . 
It is applied to a variety of areas in physics, such as cosmol- 
ogy, ultrarelativistic heavy ion collisions, or condensed matter 
physics. In describing the reheating processes during the early 
stage of inflation, it becomes necessary to trace the time evo- 
lution of slowly evolving inflaton field and quantum fluctua- 
tions ifTTUTill . In condensed matter physics it can be applied 
to Bose-Einstein Condensate (BEC), since 2PI approach is a 
candidate with properties of gapless excitation and conserva- 
tion laws ifisll . The 2PI approach with the KB equation would 
be useful also in understanding the early thermalization pro- 
cesses towards quark-gluon plasma formation in high-energy 
heavy-ion collisions [la], while it may be necessary to com- 
bine the classical field dynamics leading to the linear rise of 
entropy from the chaotic nature of the system iHtIi . 

In this paper we focus on the 2PI approach to the scalar 
0{N) theory, where identical bosons interact one another in 
the symmetric phase. The 0{N) model has been employed in 
the vacuum and thermal equilibrium analyses and it is also 
applied to time dependent phenomena j^. It can be used 
in inflationary processes of the early universe, the formation 
of Bose-Einstein condensate in the laboratory, and the chi- 



ral phase transition in heavy ion collisions. The dynamics of 
the chiral phase transition following the expansion of a quark- 
gluon plasma produced in relativistic heavy ion collisions has 
been analyzed byan 0(4) a model in the leading order of 
1/iV expansion lll8l[l9ll . 

In the analysis of the 0{N) model, the large approxi- 
mation (or the 1/A^ expansion) has been applied since 1960s 
or 1970s in both statistical mechanics and quantum field the- 
ory ||20| - |22|] . The 1 /N expansion has the advantage over the 
loop expansion that it is not restricted to small couplings. 
However a naive diagrammatic 1/A^ expansion breaks down 
when it is applied to the dynamical simulations, because a sec- 
ular (unbounded) time evolution prevents the description of 
the late-time behavior of quantum field. As an improved treat- 
ment, we can take account of the self-consistency in the 1 par- 
ticle irreducible (IPI) effective action, which corresponds to 
the first Legendre transformation of the generating functional 
with respect to the expectation value of the field, but the IPI 
effective action is also plagued by the secular problems and 
unitarity violation; for example {(j){x)'^) can become negative 
in the late time behavior ll23il24ll . Therefore in order to extract 
the stable evolution without the above secularity we need to go 
beyond the IPI technique, so that we adopt 2PI technique as 
the simplest example. Recently a systematic expansion 
of the 2PI effective action has been applied to a scalar 0{N) 
in the symmetric and broken phases 12511. By resuming all 2 
particle reducible diagrams with higher order secularity with 
respect to time, 2PI approach with the KB equation realizes 
the controlled time evolution without the secular problem and 
the unitarity violation in far-from-equilibrium dynamics lf26ll . 

Starting from controlled numerical analyses of the next-to- 
leading order (NLO) skeletons of 1 /N expansion of the 0{N) 
model in 1 + 1 dimensions with vanishing classical field in 
Ref. 1I26I1 . simulations in the broken phase in Refs. lfl2l [Till 
and estimation of nonthermal fixed point |27t| have been also 
done. It is numerically found that we need to include the 
NLO contribution of the expansion in order to describe 
quantum scattering and thermalization. Recently NNLO of 



2 



the 1 /N expansion is derived psl] and the rapid convergence 
property for moderate values of is observed ll29ll . Nu- 
merical analyses show the thermalization of the distribution 
function derived from the two-point Green's function for the 
truncated self-energy. However there is no consideration of 
kinetic entropy and its H-theorem based on the KB equa- 
tion. Thus we concentrate to derive the analytic expression 
of the entropy current in nonequilibrium from KB equa- 
tion and to show that the H-theorem is satisfied at the level of 
the Green's function analytically for the NLO self-energy of 
the expansion in the symmetric phase. In the end, we aim 
to give a criteria whether thermalization occurs or not in off- 
shell propagation for given self-energy or collision processes. 
In (j)^ theory with the NLO self-energy of the coupling ex- 
pansion, numerical simulations without gradient expansion in 
Refs. ll30l - l34ll shows the resultant thermalization. The proof 
of the H-theorem in Refs. IbsI - IbtIi is consistent with these nu- 
merical results. The work given in this paper is an extension 
of our recent work on the H-theorem in (fA theory to 0{N) 
theory IstIi . We show the H-theorem for local parts of NLO 
self-energy of the 1/iV expansion in the 0{N) theory. We 
find that the resultant H-theorem is consistent with numerical 
analyses performed for the 0{N) theory. Judging from this 
consistency, the proof of the H-theorem might be a criterion 
to confirm whether thermalization occurs or not without nu- 
merical simulation. It might be possible to investigate qualita- 
tive properties of entropy production in quantum field theories 
such as non-Abelian gauge theory, which we still have diffi- 
culty in numerical estimate, for example, due to singularities 
in massless quantum fluctuations. 



We also give an analytic expression of the entropy density 
S in thermal equilibrium derived from the thermodynamic po- 
tential (2PI effective action) with NLO skeletons of the 1/iV 
expansion. The present entropy density s° based on the H- 
theorem should correspond to the entropy density S in ther- 
mal equilibrium. We find that entropy density s'^ is consistent 
with the entropy density S in thermal equilibrium if we can 
regularize skeletons in higher loop order If singularities re- 
main, memory correction terms as shown in Ref. ifssll might 
appear. Hence we have to estimate the singularities of skel- 
tons carefully order by order in practical appUcations. 



In this paper first we review the Kadanoff-Baym equation 
with the NLO self-energy of the 1/N expansion. Next we 
derive an analytic expresson of the entropy current for the rel- 
ativistic field theory in the symmetric phase, and show the H- 
theorem for the 0{N) theory in NLO of the expansion. 
In the derivation we adopt the gradient expansion with respect 
to coordinate space, which is necessary to violate time rever- 
sal invariance of the KB equation. When we ignore higher or- 
der derivative terms, spatial fluctuations in a small volume are 
assumed to vanish. Thus the gradient expansion corresponds 
to a coarse graining procedure. In the end we compare our 
entropy density for nonequilibrium systems with entropy 
density S in thermal equilibrium. 
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FIG. 1: Integral path in the time coordinate. 



II. KADANOFF-BAYM EQUATION WITH NLO 
SELF-ENERGY OF 1/iV EXPANSION 

In this section, we briefly review the Kadanoff-Baym equa- 
tion with the NLO self-energy of the expansion in the 
0{N) model. We follow the notation in Ref. Let us con- 
sider a relativistic real scalar field 0^ (a = 1, • • -N) with the 
0{N) symmetric action. 



SI 
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X 
4L/V 



(1) 



where summation over indices is implied. We adopt a closed 
time path C along the real time axis, the path from to to oo and 
from cxD to to as shown in Fig.[T] in order to trace nonequilib- 
rium dynamics The 2PI effective action F with vanish- 
ing mean field {(/)) = (symmetric phase) is written as 



r[G] = ^Trln(G) ' + '-TrG^^G + \^G] 
with the full two-point Green's function 



(2) 



Gab{x,y) = {Tc(l)a{x)(j)b{y)) 

= 5abG{x,y). (3) 

Here Gab{x, y) defined on the closed time path C can be also 
written in a matrix notation 



Gab{x,y) 



GlUx,y) Gll{x,y) 
GlUx,y) Gll{x,y) 



(4) 



with 



Gll{x,y)={Mx)My)), (5a) 

Gilix,y)={My)Mx)), (5b) 

Giiix, y) =9ix^ - y')Gli{x, y) + - x^)Gll{x, y), 

(5c) 

Gllix, y) ^9iy' x')Gll{x, y) + 9{x'' - y')Gllix, y). 

(5d) 

The upper label "1" represents the path from to to oo and "2" 
represents that from cxd to to in the closed time path contour 
C. The term $[G]/2 contains infinite series of 2PI diagrams 
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FIG. 2: NLO skeleton diagrams of expansion in symmetric 
phase lH. 



whose lines are given by the full Green's function G. The 
stationary condition for the effective action 



5G 



= 



(6) 



gives rise to the SD equation for the non-equilibrium Green's 
function G{x, y), i.e. the KB equation, 



G-\x,y) = G^\x,y)-Y.{x,y) 



(7) 



The proper self-energy (S) and the free Green's function (Go) 
are defined as 



S =i6^G]/5G , (8) 
iGo\x,y)=-{dl+m')Sc{x-y) , (9) 

where Sab in front of G and E can be factorized. 

The NLO diagrams of expansion in $[G]/2 is given 
by the series of diagrams shown in Fig|2l which can be written 
in equations as I|25lI26Ii . 



l-flG] = J^Gaa{x,x)Gbi>{x,x) + ^Tr In i3[G] , 

(10) 



B{x,y]G) = 5l^^^\x -y) + ^ Gab{x,y)Gab{x,y) 



(11) 



The two-point correlation function G{x, y) in the 
Schwinger-Keldysh formalism is expressed by two indepen- 
dent functions; the spectral function p{x, y) and the statistical 
function F{x,y), 



G'\x,y) =Fix,y) + -p{x,y) 
G^\x,y) =F{x,y) ~ ^p{x,y) 



(12) 



(13) 



Notice that F{x,y) = F{y,x) and p{x,y) ~ —p{y,x). 
The spectral function p{x,y) contains the information about 
which states are the most realized, while the statistical func- 
tion F{x, y) determines how much a state is occupied. Then 
KB equation (|7]i can be reexpressed by using F{x, y) and 



p{x,y) as 



dr. + m'^ + X ^„t.^ F{x, x) 



6N 



F{x,y) 



tpix, z)F{z, y) - tpix, z)p{z, y) 



dx + 171^ + A 



N + 2 
6N 



F{x,x) 



P{x,y) 



dz d z'Ep{x, z)p(z,y) 



(14) 



(15) 



Here the self-energies are also decomposed into the statistical 
and spectral parts. 



^^^{x, y) =^f{x, y) + -^p{x, y) 
^'^^{x,y) ^tpix^y) - ^J:p{x,y) 



(16) 
(17) 



These are given by 



^Fix,y) =F{x,y)DFix,y) - -p{x,y)Dp{x,y) , (18) 
^p{x,y) =p{x,y)DF{x,y) + F{x,y)Dp{x,y) . (19) 
The functions Dp and Dp are given by resummation: 

DFix,y) = -IlF{x,y) 



dz 



Dpix,y) 



^d'-z (Df{x, z)flp[z, y) ~ Dp{x, z)flF(z, y) 



3N 



(20) 

-np{x,y)+ / dz°d'^zUp{x,z)Dp{z,y) 

Jto 

(21)' 



with the functions tip and lip given as. 



lipix^y) 



A 



A 



F {x,y) - {x,y) 



^pix,y) =-F{x,y)p{x,y) 



(22) 
(23) 



Here notice that 11 p and lip represent the chain part of the dia- 
grams in Fig. |2] In the limit of iV -> oo the R.H.S. of Eqs. ( fT4l ) 
and ( fTSl ) vanish and neither thermalization nor damping of 
unequal-time Green's functions occur. 



III. H-THEOREM FOR 0{N) SCALAR FIELD THEORY 
WITH NLO SELF-ENERGY 

A. Introduction of kinetic entropy 

In the derivation of kinetic entropy from the Kadanoff- 
Baym (or equivalently Schwinger-Dyson) equation (|7]), we 
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FIG. 3: Local and non-local part of the self-energy. 



adopt the G^^ and G^^ as independent functions in the 
Green's function. We start from the Kadanoff-Baym equa- 
tion (|2ll. Multiplying G from the right and left hand sides 
of Eq. (|2l), respectively, and by using G^^(a;, z)G{z, y) = 
S{x — y), we obtain, 



A(jV + 2) 
6A^ 



G''\x^y) 



i / dzi:Zny{^, z)c^''G''\z, y) = ic'^'Six ~ y) , (24) 



'd'y+rn^ + ^^^G^\y,y) 



G''\x,y) 



i J dzG'^'^ix, z)c^''^tA^, y) = zc^H{x - y) 



(25) 



The integral over z^ is from to oo, and the sign factor c"'' = 
diag(l, —1) reflects the direction of integral on branches 1 and 
2. The local part of the self-energy, shown in the first graph in 
Fig. |3] is included in the first term as the modification of the 



mass, and the non-local part Enoni is defined as the rest. 



Y.'^\x,y) 



-G''\x,y)dix-y)5ab 



(26) 



In deriving the H-theorem, we invoke the gradient expan- 
sion in the the "center-of-mass" coordinate X — {x + y)/2, 
while the relative coordinate x — y for two variables x and y 
are Wigner (or Fourier) transformed. Namely, the X depen- 
dence of the Wigner transform, 

A^{X,p)= j dse'V'A{X + s/2,X~s/2), (27) 

is assumed to be small. The gradient expansion of the Wigner 
transform of the product is given by using the Poisson bracket 
as. 



[AB]^{X,p) = J dse^' j dzA{X + ^,z)B{z,X 



A^{X,p)BUX,p) + - {A^,B.^}p^ + 0{-^) 



dA^ dB^ dAyj dB^ 



dpi^ dX^ dX^ dpt^ 
(Here after we omit the subscript w.) 



(28) 
(29) 



We take the difference of Eqs. (l24l) and dZSl l and carry out the Wigner transformation. By assuming that to is small enough 
(tp — cx)) or X'^ is large enough, we find 



2ip- 



d i X{N + 2) [ d<^+^k (dG'^'iX^k) dG^\X,k)\ d 

dX ) ' 



G-\X,p)^iI-\X,p) 



dX 2 6N J (27r)''+i V dX 

I^\X,p) = [S^UC^'G'*' - G'^^C^'^SnonlG'^'']^ 

= j dse'P-' J dz + s/2, zyG'^'iz, X - s/2) - G^^X + s/2, z)c='^Ef,„i(z, X - s/2)] 



(30) 



(31) 



where p is the momentum conjugate to x — y. The gradient expansion with respect to the center-of-mass coordinate X is 
appropriate when the X-dependence of the system is smooth enough (See for example ifssl [39|] ). The 1st order terms in the 
gradient expansion are taken for the Green's functions and the self-energies. We assume here that the interval of a;" — y° can 
take sufficiently large values for a fixed X^ . 

By using Eq. ( |28] ). the gradient expansion of /"'' for (a, 6) = (1, 2) and (2, 1) is found to be. 



^nonl'"-^ JPB I '^nonl/pB 



= - C{X,p) + {Re^n, G^^}^^ - {RcGr, Si^mjpB 
=- G(X,p) + {ReSfl,G2i}p3 - {RcGfl,E2i^,i}p3 

I 



(32) 
(33) 



where the term G, 

CiX,p) = t F{X,p) - ^f{X,p) p{X,p)] (34) 

may be identified as the collision term in the Boltzmann limit. 
In the last line of Eq. ( |32] |. we have utilized the following 



relations, 

^12,21 



F±-p. 
2^ 



gii,22 -^^^^ ^ 5.11^22 ^^^^ -j^ej.^ 



(35) 
(36) 
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The equation of motion for G" (a = 12, 21) are then rewritten 



as 



{h, G^^ + {RcGij, S^onilpB 



V21 
^non 

2 A(iV 



h{X,p) =p^ - 



= -c, 



(37) 
(38) 



6iV 



ReSij(X,p) 



(39) 



Here we notice that the leading order of the gradient expan- 
sion of G(X,p) is already 1st order - 0{d/dX) Sill. 

Next let us take the difference of Eq. dJTb multiplied by 
In(iGiVp) and Eq. (O multiplied by \n{iG'^^ / p). Then 
we integrate the resultant expression over (f^+^p/ (27r)'^+^ to 
arrive at the following equation: 



(40) 



In Eq. (|40| ). we have defined the entropy current s^(X) as 



d+l 



(2^) 



1 aRcs^ 

2 



-G^^ln. 
d_ 

d 

dptj. 



Ep iG 



In 



P 

iG2i 



(41) 



where we have used the relations i(S]^^ — E^^) = 2RcE/j and 
z(G" - G22) = 2RcGfl fromEq. We have also adopted 
the approximations, 



e.g^Vp, e^i^e.g^Vp^ 



(42) 



in the 1st order gradient expansion 11351 136ll38ll . 

In the end we shall write the two-point functions in the 
form of the Kadanoff-Baym Ansatz G^^ = —ipf and G^^ = 



J 



—ip{l + f) with a real function /, then the above entropy 
current can be rewritten as 



P 



(27r) 



d+l 



p / 1 aRcE 
i V 2 dpf. 



Ep 1 dReGn 



i 2 dp^ 



where we introduced the notation 

a(X,p) = -/ln/ + (l + /)hi(l + /) 



(43) 



(44) 



In the quasiparticle limit Enoni ^ 0, we know that the Green's 
function is given by the occupation number rip, 

Gi2 = _ ipf = 2ttS{{p")^ - nlm-p'') + Tip) , (45) 

G^i ^ - tp{l + /) = 27r<5((/)2 - nl){9{p°) + np) , 

(46) 



X{N + 2) f dfik 
ON J {2ttY 



(47) 



Then we reproduce the normal entropy current for bosons 

d'^p pt" 



[-Up hi Tip + (1 + Tip) ln(l 



(48) 



The remaining work is to prove that the R.H.S. of Eq. ( |40l ) 
becomes positive definite. 



B. NLO self -energy of 1 /N expansion and H-theorem 

In this subsection we prove the H-theorem for the NLO self- 
energy in the 0{N) model by use of the gradient expansion 
of the Kadanoff-Baym equation. In the previous section and 
Ref. fyh . we find the relation ( l40t where the divergence of 
entropy density is proportional to G = i(EpF — Ei?p) for the 
scalar field theory. We shall Wigner (or Fourier) transform 
Ei? and Ep in addition to Dp and -Dp with 11^ and tip. 



First in order to derive Dp{X,p) we calculate Dii{X,p) and Da{X,p) by multiplying 6{x'^ 



yO) and -9{y° - x°) in (EB 



6 



and Wigner transform. Then Dfj {X, p) can be calculated as. 



DRix,y) =6{x'' ^y'')bp{x,y) 
X 



— - y°)n,(x, y) + —e{x' -y°) I dz Xl^ix, z)Dp{z, y) 

^e{x° y°)U,{x, y) + ^0(^0 _ ^o^ y np(x, z)Dp{z, y){9{x'> z") - 0(y" - z")) 



(due to 6l(x" - y")9{x° - z") = e{y" - z")6'(a;° - y") + Oi^' - y°)e{x^> - z")) 
_ —0(j:O _ yO)n y) + — 



to 



11b X, w H 

3iV '^-^ 3iV 



dz ri,(x, z)^p(z, y)0{z^ - y")0(a;O - z") 
dz Unix, z)Dii{z,y) . 

I 



(49) 



By taking to — s> — oo, Wigner transformation, and using the a similar way Da{X,p) = — F.T.6'(y" — x^)Dp{x, y) can be 
gradient expansion relation ( |28] |. Dii{X,p) can be written in written as 
the leading order of the gradient expansion as 

A 



Dr{x,p) = -^{nR{x,p) - tIr{x,p)Dr{x,p)) .(50) 



bA{x,p) = - A \tiA{x,p) - \iA{x,p)bA{x,p) 



(52) 



As a result we obtain Dji{X, p) can be derived as 

Dr{X,p) 



^^i\R[x,p) 
1- 3^nfl(x,p) 



where we have used the relation 8{y^ — x^)0{y^ — z*^) = 
e{x° - z°)9{y° - x°) + 9{z° - x°)e{y° - z°). As a resuh we 
obtain 



(51) 



DAiX^p)^^^"^ 



where we have taken only the leading order (local parts) of the 
gradient expansion in the convolution by use of ll45ll In 



J 



1 - ^IlAiX,p) 



(53) 



By using the relations bp{X,p) = bii{X,p) — Da{X,p), ( BTl ) and ( [53] l. we obtain 

bpix,p) ^bR{x,p)~bA{x,p) 



-^%{X,p) 



where we have defined the effective coupling as, 

\cs{X,p) = 



(l-^t\R{X,p)) (l-^t\A{X,p) 



A 



1-^\Ir{X,p)] l_^fl^(X,p) 



^^np(^,.), (54) 



(55) 



In a similar way Df{x, y) can be derived by rewriting as 



bpix^y) 



t\p(x, y) 



dz [bpix., z)np(z, y) - bpix, z)IIf{z, y) 

A - A l"^ / - 

= - ■^ni?(a;,y) + — j dz ybF{.x,z)\lA{z,y) + bR{x,z)ilF{z,y) 

Then the Wigner transform Df{X, p) is derived by taking — oo as 

Ks{x,p) 



(56) 



bF{x,p) = — —ivf + — 



DfUa + DrUf 



(l - ^Ur) (l - 3^11^ 



3A^ 



IlFiX,p) , (57) 
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where we have used Eq. (ISTT i. By use of ( l54b and ( fSTl ) we obtain the following relations with respect to Wigner transformed 

tF{X,p) and tp{X,p) from ([Till and Gil as 



dq { F{X, q)DF{X,p- q) - -p{X, q)Dp{X,p- q) 

_Acff(|^^^\ (p^x,q)iiF{X,p-q)-\p{X,q)iip{X,p-q) 



^FiX,p) = 



tp{X,p) = I dq i^piX, q)DF{X,p- q) + F{X, q)Dp{X,p - q) 

dq (-^^^^l^f^) q)iiF{X,p - g) + F{X,q)il,{X,p - q)) , 



(58) 



(59) 



where IIf{X,p) and Ilp{X,p) can be expressed from (l22l i and (|23]) as 



Uf{X,p) =^ y (^F{X,p - q)F{X, q) - -^p{X,p - q)p{X, q) 
ilpiX.p)^^ 1 dqF{X,p-q)p{X,q) . 



(60) 
(61) 



In the end T,f{X, p) and f^p{X, p) can be reexpressed by use of ( l60b and dMT i as 



( g)^F(^,P - 9) - q)Dp{X,p~ q) 
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Acff(X,p - g) 
3iV 



F{X,q)F{X,l)F{X,p-q-l) 



-^F{X, q)p{X, l)p{X, p-q-l)- q)F{X, l)p{X, p - q - I) 

tp{X,p)= / dqdl \p{X,q)bF{X,p - q) + F{X,q)bp{X,p - q) 



dqdl 



Xcs{X,p - q) 
3N 



p{X,q)F{X,l)F{X,p-q-l) 



-p{X, q)p{X, l)p{X, p-q-l) + 2FiX, q)F{X, l)p{X, p-q-l) 



(62) 



(63) 



Finally we obtain the following colHsion term C{X,p) = i{'Ep{X,p)F{X,p) - T.F{X,p)p{X,p)) in Eq. ( |34l ). 



187V 



J dp2dp3dp4 S'^'^^ipi +P2-P3 - P4)Ks{X,p+ P2)K{X,Pi,P2,P3,Pa) , 



(64) 



KiX,pi,p2,P3,Pi) =G^\X,p,)G^^iX,p2)G^\X,p3)G^\X,p4) - G^\X,pi)G^\X,p2)G^\X,p3)G^^iX,p4) , (65) 
where we have used the self-energy (|62] |. ( l63T l, — g) = g), — = — /'(-'f, q) and 



p{X,l)F{X,r)+F{X,l)p{X,r) = i [G^\X , l)G^\X , r) - G^^ {X , l)G'^ {X , r)] , 



(66) 



with the relations = i [G2i(X,p) + Gi2(X,p)] andp(X,p) =i [G^\X,p) -G^^{X,p)]. 
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As a result, the divergence of the entropy current per particle component can be written as 



d^s^iX) = - / dp In 



A 



36N 

A 
367V 

A 

'367V 
A 

'367V 
X hi 



tpiX,p)F{X,p) - tF{X,p)p{X,p) 

G''{X,p,) 
G^HX,p,) 

G^\X,p,)G^\X,p^) 



G^^iX.p) 
J dT K{X,pi,p2,P3,Pi) In 

• ^ J dr K{X,pi,p2,P3,p4) In 

■'-JdrKiX,p„p2,p„p,) In Lg-(X,pOG-(X,p,)G-(X,P3)G-(X,P4) 
■\Jdr [G^\X,pi)G^^X,p2)G'\X,p3)G^\X,p4)~G^\X,pi)G^\X,p2)G^^{X,p3)G^^{X,pi)] 

f G^^{X,p,)G^^X,p2)G^HX,p3)G^HX,p,) 
\ (X, p^)G^^ {X, P2)G12 (X, P3)G^^ {X, P4) 



G21(X,pi)G21(X,p2 

G^'{X,p,)G^HX,p2)G'HX,P3)G'HX,p,) 



> 



dr = dp.,5'^+^{pi +P2- P3 - Pa)Ks{p + q) ■ 



(67) 
(68) 



In this derivation, we have utilized the inequality (x — 
y) \n{x/y) > 0, and the symmetry of K under the exchange 
of momenta, pi <-!> p2 and {pi,p2) ^ {P3,P4), 

K{X,pi,p2,P3,P4) =K{X,p2,Pl,P3,P4) 

= - K{X,p3,p4,pi,p2) ■ (69) 

Hence we have proved the H-theorem for the NLO of ex- 
pansion in the range of 1st order gradient expansion. We can 
see that any change of Green's functions with nonzero colli- 
sion term G = i SpF — J^pp 7^ contributes to entropy 
production. 

The equality in Eq. ( |67| i holds when K = is satisfied 
for all combinations of momenta which satisfy the energy- 
momentum conservation. This condition is realized when 
G^^ / G^2 is a linear function of p. 



f{X,p) 



G^^{X,p) 



fiX,p) 



= a{X) + P''{X)p^ (70) 



where a{X) and I3{X) are arbitrary functions for the center of 
coordinate X. Thus we notice that at equilibrium the famous 
local equilibrium distribution function 



fiX,p) 



5a(X)+/3M(X)p^ _ I 



(71) 



is reaUzed. This is one of our new results. 

Equation i64[ shows that C actually has the form of the 
collision term of bosons. The first term in K represents the 
loss term, which is proportional to the product of probability 
G^2(pi)G^2(p2) oc f{pi)f{p2) for particles having pi and 
P2, and proportional to the bosonic enhancement factor in the 
final state, G-'Hp3)G^Hpi) « (1 + /(p3))(l + fiPi)), as 
schematically shown in Fig. ID 



G^\X,l) G^\X,r) 

\/ 

lD^.s{X,p + q) 

G^'^{X,p) G^\X,q) 
FIG. 4; Collision term C 



IV. ENTROPY FROM THERMODYNAMIC POTENTIAL 
IN THERMAL EQUILIBRIUM 

In this section we derive an analytic expression of the en- 
tropy density in thermal equilibrium from the thermodynamic 
potential. Its derivation is based on Refs. ll40l l4lll and we 
adopt the imaginary time formalism. The new point in this 
article is the derivation of entropy density in the presence of 
NLO skeletons of the expansion at equilibrium. At ther- 
mal equilibrium, the entropy density is given by differentiat- 
ing the 2PI effective action or thermodynamic potential with 
respect to temperature T. We find that the resultant entropy in 
the 0{N) theory corresponds to the kinetic entropy ( |43] | with- 
out memory corrections in the limit of thermal equilibrium for 
NLO skeletons of the 1 /N expansion. 

First let us begin with the Lagrangian 



4!iV 



(72) 



where a = 1 A^. The thermodynamic potential = —PV, 
which corresponds to the 2PI effective action, of the scalar 
field can be written as the following functional of the full 
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Green's function GlUH]: 

pn[G]/N = - llogZ 

= iTrlogG-i - ixrSG + \<S>[G]/N , (73) 

where Z represents the generating functional, Tr denotes the 
trace in configuration space, /3 = l/T, and $[G]/2 is the 
infinite series of 2PI skeletons. Here we have expressed the 
functional per particle components N. The self-energy S is 
defined by 



SG 



(74) 



The self-energy is related to the Green's function G by the 
Schwinger-Dyson equation G^^ = Gq^ + E, where Gq^ is 
the bare Green's function. The Schwinger-Dyson equation is 
obtained by imposing the stationary condition of the thermo- 
dynamic potential il: 



Sn[G]/SG = 0. 



(75) 



We discuss thermodynamic potential only under the stationary 
condition in the following analyses. 

The thermodynamic potential ^1[G] can be given in the mo- 
mentum space as 



n 

NV 



(76) 



where d + 1 is the space-time dimension and = 
l/{e^'^" — 1). The summation with Matsubara frequencies 
(Matsubara sum) can be converted to the integral over fc" of 
the integrand multiplied by the Bose function f{k'^) by using 
standard contour integration techniques. The Green's function 
G(fc) is written analytically in terms of the spectral function; 



G(/fcO,k) 



dq' p(g°,k) 
(27r) - fco 



(77) 



For real variable fc", we define the regularization of the inte- 
gral as. 



ImG(/c", k) =ImG(fc° + ie, k) 



p(fcO,k) 



ReG(fc°, k) =ReG(fc" + ie, k) = ReGfl(fc°, k) 



(78) 
(79) 



The regularization for the self-energy E is defined in a similar 
way. 

Next we shall derive the entropy density S by differentiat- 
ing n by T: 



S = - 



dn/v 

dT ■ 



(80) 



The thermodynamic potential, as given in Eq. (|76] l, depends 
on its temperature through the distribution function /(fc") and 
the spectral function p (or the Green's function G). By use 



of Eq. ( fTsT l, temperature derivative of the Green's function or 
spectral function p cancels out in the entropy density due to 
SD equation, and we obtain ll40l - l43ll 



S 

with 

S' = 



N 

■N 



(27r)'*+i dT 
d'^+^k df[k°) 
{2ttY+^ dT 



Im logG-i(fc) 
ImI](fc)ReG(fc) + 5' (81) 



d{T^/V) 



dT 



■N 



d'^+^k dfik°) 
(27r)'^+i dT 



RcI](fc)ImG(fc). 
(82) 



By using the relation 



dT 

a{k^) 



log 



1 + ./ 



d(j{k°) 



dk° ^ f 9fco 
/log/ + (l + /)log(l + /), 



(83) 
(84) 



we can partially integrate with respect to and obtain 



S-S' 



=N 

=N 



d'^+^k 

d'^+^k 

d'^+^k 
(27r)'*+i 



-Im G 



dG- 



dk^ 



— (ImERcG) 



^,0 9ReE\ ^ ^ dRcG 
2k" jTTT^ 1 ImG + ImS 



where we have used G 



9fco ; 

1 dRcT^R 

2 dk" 



Gn 



dk" 
I dRcGR 
2 dk" ■ 



a 
(85) 

S and 



the definitions ( iTSl l and ( |79] l. The remaining part is S' in the 
above expression. 

In the end we shall examine whether 



5' = 



(86) 



is satisfied for both LO and NLO skeletons of the 1/iV expan- 
sion. 

It is easy to verify S' = for LO skeleton of the 1/iV 
expansion. The LO skeleton is given by the expression 



91/ 94 ^ 



d'^k 
{2ny 



Giio,k) 



2V 

XN r d'^+^k, [d^+^ ^faOw.,o^ 



24 



where we have used Eq. (fTTT i and the relation 

1 



(87) 



(88) 
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Then we take account of the temperature dependence of the 
occupation number function /. By differentiating by T we 
obtain the first component in S' (|82] | for the LO skeleton as 



d f T 



dT \ 2V 



-$(LO) 



XN 



G 



(2^) 



(89) 



The second component of S' ( [82l l has the real part of the self- 
energy which is given by opening one line of the skeleton. 
The real part of the self-energy S 



A 



SG 



is written as 



-/(fcO)p(fc)(90) 



'6 J (27^)''^ 

This relation ( |90l l gives the following second component of iS' 



12 J (27r)''+i J {2t:Y+ 



where we have used Im G = p/2. We can confirm that two 
components (|9T] i and ( [89] l cancel each other precisely and that 
5' = is verified for LO skeleton of 1 /N expansion. 

Furthermore we shall derive the expression of S' for the 
NLO skeletons of the expansion in the 0{N) theory by 
writing down two components of S' explicitly and investigate 
whether 5' = is satisfied or not in this order We find the n+ 
1-loop {n > 2) order NLO skeleton per particle component is 
given by 

_ ^$("+1) 
2V 



{2ttY 



(27r)° 



2n\ 6 J ^„ . , , 
X G(w,q)G(a;i,ki) 

X G{lo + - a;2,q + ki - k2)G(a;2,k2) 
X • • • 

X G{uj + uji -a;„,q + ki - kn)G(a;„, k„). (92) 

For the above skeletons the explicit form of first and second 
component in S' is given in Appendix A. Then we find that it 
might be possible to prove S' ^ for each loop order of NLO 
skeletons of the 1 /N expansion. However we also notice that 
these skeletons might have singularities at the points where 
energy denominator vanishes for 7i > 3. Memory correction 
terms can appear when momentum integration with these de- 
nominators can not be regularized ll35il44ll . Then we have to 
carefully estimate the singularities and add the memory cor- 
rection terms in the definition of entropy current. Hence the 



possibility of the appearance of memory correction terms is 
not completely denied. 

V. SUMMARY 

In this article, we have discussed the entropy density and 
the H-theorem in the Kadanoff-Baym equation for the 0{N) 
theory with the next-to-leading order (NLO) self-energy of 
the expansion. While many numerical analyses show 
thermalization in the Kadanoff-Baym equation takes place, it 
is desired to confirm thermalization analytically for a given 
form of the self-energy. Thus we have derived the relativis- 
tic kinetic entropy and confirmed that the H-theorem is satis- 
fied for local parts of the NLO self-energy. We have adopted 
the gradient expansion and taken account of the 1st order 
terms of the Kadanoff-Baym equation. The gradient expan- 
sion is an appropriate approximation in treating moderately 
changing systems in space-time. Since we ignore higher or- 
der derivative terms, spatial fluctuations in a small volume are 
assumed to vanish. Thus the gradient expansion corresponds 
to a coarse graining procedure. Time reversal invariance is 
violated by the gradient expansion of the (time reversible) 
Kadanoff-Baym equation. Then it is suggested that the KB 
dynamics obtains the time arrow at the level of the Green's 
function, and that the observed system at late time will reach 
the maximum entropy state where local equilibrium is real- 
ized. The reliability of the gradient expansion and the form of 
the self-energy are important in the proof of the H-theorem. 
From the proof of the H-theorem we find that any change of 
the Green's function with nonzero collision term contributes 
to entropy production. The entropy production is ensured at 
the level of the Green's function (No need of quasiparticle 
approximation). When local thermal equilibrium is realized, 
the collision term contribution vanishes and entropy ceases to 
increase. The proof of H-theorem is consistent with the nu- 
merical analyses of the 0{N) model. Thus the proof might 
provide a criterion whether thermalization occurs or not even 
in the cases when numerical analyses are difficult to perform. 

We have also considered the entropy density S derived from 
thermodynamic potential, and compared with the kinetic en- 
tropy density. When we only take account of the principal 
value part in the momentum integral for S, S corresponds to 
the kinetic entropy density in thermal equilibrium. However if 
singularities remain in the energy denominators, memory cor- 
rection terms might appear Hence entropy current have the 
possibility to have additional terms in its definition. 
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Appendix A: The n-loop order skeleton and S' 

In this appendix we shall confirm iS' = for the NLO skeletons order by order in the coupling expansion. The n + I loop 
order component of $ in 1/A^ expansion is given in Eq. (|92b . This equation is transformed as 

with Eq. dTSl l. We shall sum up (j = 1, • • •, n) by use of the relations 

" . 1 1 " f \ 1 \ 1 



n 

n ,o^fcO_,^_, (/fe") + /bl-) + l) (A2) 



i=2 



with IMll and -/{-p^) = /(p?) + 1, and 



1 " 1 1 " 

" i.o _„„ IT- 



1 



" 1 1 " 1 



Then we obtain the factor 5'^ in < IAlb as 



Sec 



fj + fcO - a; - fc? + kO - ipO + fco - c.) J^^ + fc? - pO fcO ^^4) 



with /(pO + k'^){.f{p°) + f{k^j) + 1) = f{p'j)f{k'j)- In the end by use of the relation 



J ^/'jV J V^-j „0 _i_ jLO _ „0 „ 1,0 11 „0 I 1,0 _ „0 _ 1,0 



with 



1 " 1 1 " 1 

S + fc? - w - fc? " 9" -CJ n pO + fcO _k°-qO 

" 1 1 " 1 

'§ + -kl-Loq^- (p« + fcO - fc?) ^. (p? + k'l fcf) - (p^" + fcO - fc?) ^^^^ 



and 



1 ^ 1 1 



(A7) 
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we obtain 



2V 



2n V 6 



(2^ j J2^P^^^P^^ 



(=2 



nw.,o. /(fcg) + /(g°) + i A /(Pl^) + /(fcn + i 



J=2 



(A8) 



(Here for ti = 2 since + /(p^) + 1 = .f{kf)f{p^){ 



1), by use of the interchange {q, ki) o {p2, ^2) the above 



terms in the bracket can be transformed as 



,.0>/(pS) + /M) + 1 . 1 



e T 



(A9) 



so that the skeleton ^^"^^ is analytic function and the fraction with energy denominator can be written by its principal value. ) 
Hence the first component of S' is given by 



d 



T 



dT \ 2V 



6 / dT 



(2^) 



d+l 



Pi 



d+l 



(27r)'^<5('*)(pj+Pj-q-ki)p(p,)p(fc,) 



n 



/(P?) + fik?) + 1 



+ fco - pO - k° 

The second component J ^^J^^!li KeT.''"'~^^\q)lmG ^-^Q^ ^ cancels the principal value of dAlOl ) where ImG = | and 



~SG. 



ab 



XV 5, 



6 J N 



d'^k 



xG{oj + 0J1 - W2,q + ki - k2)G'(a;2,k2) • • ■ G{oj + oji - u;„,q + ki - kn)G(a;„, kn) 



-2(-^V'^ 



n 

U=2 



N J (2n)° 

d'^+^k, r d<^+^ 



{2n 



\d+l 



{2n 



-g^(2^)'^<5('^)(pj + kj - q - ki)p(p,)p(fc,) I X 



(AlO) 



(All) 



with changing variables uj — i2TrlT — > g". As a result we find iS' = for its principle values due to the cancellation order by 
order. When singularities of skeletons for n > 3 remain, memory correction terms appear in thermal equilibrium in the final 
expression of entropy density S as Refs. Issiliill . 
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